Walls inhibit chaotic mixing 
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We report on experiments of chaotic mixing in a closed vessel, in which a highly viscous fluid is 
stirred by a moving rod. We analyze quantitatively how the concentration field of a low-diffusivity 
dye relaxes towards homogeneity, and observe a slow algebraic decay of the inhomogeneity, at 
odds with the exponential decay predicted by most previous studies. Visual observations reveal 
the dominant role of the vessel wall, which strongly influences the concentration field in the entire 
domain and causes the anomalous scaling. A simplified f-D model supports our experimental 
results. Quantitative analysis of the concentration pattern leads to scalings for the distributions 
and the variance of the concentration field consistent with experimental and numerical results. 



Low- Reynolds-number fluid mixing has a variety of ap- 
plications ranging from geophysics to industrial mixing 
devices. While turbulent flows lead to highly efficient 
mixing, simple laminar flows with chaotic Lagrangian 
dynamics also promote rapid homogenization [1]. Dy- 
namical systems approaches based on flow kinematics 
have provided a first insight into chaotic mixing 0, 
A deeper understanding of homogenization processes is 
gained by examining the interplay between chaotic stir- 
ring and diffusion. Several experimental 0, IH and nu- 
merical d, 0] studies obtained an exponential decay for 
the variance of a diffusive scalar concentration field in a 
chaotic mixer. This behavior is attributed to an asymp- 
totic spatial structure of the scalar dubbed strange eigen- 
mode [8[ , that results in a global exponential decay of the 
spatial contrast. However, these theories focus on ideal 
mixing systems, e.g. with periodic or slip boundary con- 
ditions, far from the reality of industrial mixin g de vices 
with solid no-slip walls. It has been suggested [9|, [lfj that 
mixing might be slower in bounded flows, but experimen- 
tal evidence is still lacking. 

In this Letter, we study experimentally dye homoge- 
nization by chaotic mixing in a 2-D closed flow. Precise 
measurements of the concentration field yield "slow" al- 
gebraic decay of an inhomogeneity, at odds with the ex- 
pected exponential decay. We relate quantitatively this 
slow mixing to the chaotic nature of trajectories initially 
close to the no-slip wall, which end up escaping in the 
bulk and slow down the whole mixing process. 

A cylindrical rod periodically driven on a figure-eight 
path gently stirs viscous sugar syrup inside a closed ves- 
sel (Fig. [1] (a)). The stirring scale is comparable to the 
vessel size, in contrast to other devices such as arrays 
of magnets 0, [f| . This protocol is a good candidate for 
efficient mixing: we can observe on a Poincare section 
[ill ] (Fig. [2 (a)) - computed numerically for the cor- 
responding Stokes flow - a large chaotic region span- 
ning the entire domain, including the vicinity of the 
wall. The signature of chaotic advection can also be ob- 
served in Fig. [T] (a), where a complex lamellar pattern 



is created by the stretching and folding of an initial dye 
blob into exponentially thin filaments. The fluid viscos- 
ity v — 5 x 10~ 4 m 2 • s _1 together with rod diameter 
I = 16 mm and stirring velocity U — 2 cm ■ s _1 yield a 
Reynolds number Re = Ul/v ~ 0.6, consistent with a 
Stokes flow regime. A spot of low-diffusivity dye (Indian 
ink diluted in sugar syrup) is injected at the surface of 
the fluid, and we follow the evolution of the dye concen- 
tration field during the mixing process (see Fig.[T]). The 
concentration field is measured through the transparent 
bottom of the vessel using a 12-bit CCD camera at reso- 
lution 2000 x 2000. 

Despite the exponential stretching occurring in the 
bulk, the resulting variance tr 2 (C) of the concentration 
field (measured in a large central rectangular region) de- 
cays surprisingly slowly with time t as a power law t~ m 
with m ~ 3.2 (Fig. Q] (b)), and not exponentially as ex- 
pected. This behavior persists until the end of the exper- 
iment (35 periods), by which time the variance has de- 
cayed by more than three orders of magnitude. Moreover, 
concentration probability distribution functions (PDFs) 
shown on Fig. [T] (c) exhibit wide power-law tails on 
both sides of the most probable value. The probability 
of "white" (zero) concentration decays very slowly with 
time, whereas the peak shifts towards lower values. 

In order to understand these surprising scalings, we 
first describe the various mechanisms at play during the 
mixing process. A blob of dye, initially released close 
to the vessel center, is transformed into a complicated 
pattern expanding towards the wall with time. We dis- 
tinguish at each instant the growing "mixed region" , de- 
limited by the advected blob frontier, and the remaining 
wall region where C = 0, in the vicinity of the vessel 
wall. This distinction is obvious in Fig. [1] (a) where one 
can observe a central heart-shaped mixed region and an 
annular unmixed wall region. As the chaotic region spans 
the entire flow, fluid particles initially close to the zero- 
velocity wall eventually escape from the wall region to 
wander through the whole chaotic region. Trajectories es- 
cape along the unstable manifold of parabolic separation 
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FIG. 1: (a) Chaotic mixing experiment in a closed vessel: 
a rod moves periodically on a figure-eight path (see Fig. [5] 
(a)) and transforms an initial spot of dye into a complicated 
filamentary pattern. Inset: Poincare section obtained numeri- 
cally for the corresponding Stokes flow. Note that the chaotic 
region spans the entire domain, (b) Evolution of the vari- 
ance of the concentration field in a fixed central region. ■: 
experiment , o: numerical simulation (see Fig. [3] for a de- 
scription). Solid line fits: contribution of the white pixels 
(7\y — (2 log t + log wb)/ (log F x t 2 ) derived below, (c) Exper- 
imental concentration PDFs in the bulk at time t = 13, 17, 31 
periods. Both sides of the peak can be fitted by power laws 
{Cmax — C)~ 2 (red and blue plots). Inset: left ("light-gray") 
tail of the peak, P(C) against \C ma x(t) — C\. 



points on the wall 0, [l2| ■ The signature of such an escape 
path can be visualized on Fig. [1] (a) where unmixed fluid 
from the wall region is "sucked" inside the mixed region 
through its white cusp, close to the rod. This results in 
the periodic injection of broad white strips that can be 
observed inside the mixed region. The mixed region then 
grows towards the wall to make up for this mass injection. 
Incompressibility combined with zero-velocity condition 
at the wall leads [9] to a shrinking distance between the 
mixed region border and the wall scaling as d{t) oc t . 
This scaling is verified experimentally. The area of un- 
mixed fluid from the wall region injected at each period 
inside the mixed region then scales as d(t) oc t~ 2 . As a 
result, the mean concentration value inside the bulk de- 
creases with time as (1 — <i(i)) _1 . Simultaneously, the 
mixed region is stretched and folded at each half-cycle 
of the rod movement (see Fig. [5] (a)) in a baker's-map- 
like fashion 
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FIG. 2: (a) Transport mechanisms: (i) the rod stretches and 
folds the mixed region and (ii) a white strip of unmixed fluid 
is injected between the two folded parts, (b) Stretched strips 
of dye (black) are smeared out by diffusion on a scale wb 
(gray). The concentration C of a pixel x is then given by 
adding the concentrations coming from strips inside a box of 
size wb around x. 

but separated by the newest injected white strip. Also 
note that the part "attached" to the rod has experienced 
more stretching than the one "left behind". Briefly, (i) 
chaotic stretching imposes that the typical width of a 
dye filament in the bulk shrinks exponentially down to 
the diffusion or measurement scale, yet (ii) wide strips of 
unmixed fluid of width d(t) oc 1/t 2 are periodically in- 
serted between these fine structures. In the following we 
derive how these two effects lead to the observed scalings. 

For this purpose, we simplify the 2-D problem by char- 
acterizing only 1-D concentration profiles C(x,t) along a 
secant to the stretching direction - the dashed segments 
on Fig. [3 (a) - that is, we neglect the variation of the con- 
centration along a dye filament on a scale comparable to 
the vessel size. We thus call from now on "mixed region" 
the intersection of the 2-D mixed region with such a seg- 
ment. The effect of the mixer during a half-period then 
amounts to the action of a one-dimensional discrete-time 
map that transforms concentration profiles by inserting 
an interval of width d t of fluid from the wall region be- 
tween two inhomogeneously compressed images of the 
mixing region at the previous time (see Fig. [2] (a)). Such 
a map /, defined on [0, 1] for simplicity, evolves concen- 
tration profiles as C(x,t + 1) = C(f~ 1 (x),t) and meets 
the following requirements: (i) it is a continuous double- 
valued function to account for the stretching/folding pro- 
cess; (ii) x — is a marginally unstable point of the 
correct dynamics close to the wall are indeed reproduced 
by imposing f~ 1 (x) ~ x + ax 2 + ■ ■ ■ , a > for small 
x; (iii) because of mass conservation, at each x, the local 
slopes of the two branches add up to 1 . Other details of 
/ are unessential for our discussion. Diffusion is mim- 
icked by letting the concentration profile diffuse between 
successive iterations of the map. This model is a mod- 
ified baker's map 
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13]. However, it should be noted that the 



two folded parts are not stacked directly onto each other 



with a parabolic point at x = 0, 
whereas the dynamics are purely hyperbolic in a classical 
baker's map. Numerical simulations for a specific choice 
of / lead to results shown in Fig. Q] (b) and Fig. [3j both 
the power-law evolution of the variance and the different 
aforementioned features of the experimental PDFs are 
reproduced by the simulations. 
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FIG. 3: Numerical concentration PDFs (t = 15, 18, 25) ob- 
tained by letting concentration profiles evolve as C(x,t + 1) — 



C(f-\x),t), f(x):f 1 (x) = x- 



+ (7-l + a)x J ;/ 2 (a;) 



1 — aar + (7 — 1 + a)x 6 with a = 0.9 and 7 = 0.55. Up- 
per (resp. lower) inset: light (resp. dark) gray tail, P(C) 
against \C ma x — C\. We observe the same power-law decay 
(Cmax ~ C)~' 2 on both sides of the maximum as in the exper- 
iment (Fig.Q](c)). 



The map transforms an initial blob of dye of width 
So into an increasing number of strips with widths 
so7i'''7tj resulting from different stretching histories 
inside the mixed region, where j t is the compression 
experienced at time t. White strips also experience 
this multiplicative stretching from their injection time. 
As the mixed region grows towards the wall, different 
values of stretching will be sampled. It will never- 
theless be justified below that the concentration mea- 
sures only require knowledge of the stretching histories 
traced back during a finite number of periods. This 
allows us to define in a "quasi-static" approximation 
a slowly varying instantaneous "Lyapunov exponent" 

T(t) = exp((- log (| df d } x) \))mr), that is the geomet- 
ric mean of the compression taken over the mixed region 
(MR) at time t. Note that as the white strips are in- 
jected close to the center of the domain, the two branches 
of / have comparable mean slopes, yielding the estimate 
T ~ 0.5. Without diffusion, dye strips would have a 
typical width sqT 1 at time t. However, the balance be- 
tween stretching and diffusion imposes that the width of a 
strip stabilizes at the Batchelor scale wb — \J — T 2 ), 
where k is the diffusion coefficient, wb is thus the small- 
est lengthscale that can be observed in the concentra- 
tion profile, and different elementary strips may overlap 
(Fig. [2] (b)). Since the concentration field is probed by 
averaging it on the pixel size w px , which is smaller than 
Wb, the concentration at a pixel is given by adding the 
contributions from strips contained in a box of size Wb 
around the pixel. Hence we characterize P(C) by con- 
sidering the different combinations of strip widths for a 
zero-diffusivity dye - the concentration profile on Fig. [2] 
(b) - that one might find in a box of size wb- We will 
distinguish between three generic cases corresponding to 



three different regions of the histogram P(C) (see Fig. Q] 
(c) and Fig. [3]): a white (W) peak at C — correspond- 
ing to injected white strips still wider than wb, light gray 
(LG) and dark gray (DG) tails corresponding to respec- 
tively smaller and larger concentrations than the peak 
(mean) concentration. Once we have quantified the pro- 
portion of boxes contributing to these different values of 
C, the variance will be readily obtained as 



a 2 (C) = {C- {C)) 2 P(C)dC = a 2 w + a\ G + 
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Let us start with white (zero) concentration measures 
that come from the stretched images of white strips in- 
jected before t. White strips injected at an early time 
have been stretched and wiped out by diffusion, that is 
their width has become smaller than wb ■ Hence the old- 
est white strips that can be observed have been injected 
at the time ti(t) such that d ti T t ~ ti — %. Note that 
from £ w hite defined by dt whitc = wb, the injected white 
strip is smaller than wb and no white pixels can be ob- 
served. Before i w hito, the number of white pixels is pro- 
portional to nw — ^ti(t) — d t ct (t — ti)(tit)~ l for large 
t (using d t cx i -1 ). As t — U ~ (21ogt + logics)/ logT, 
nw — (2 log* + logWB)/(logr x t 2 ). We deduce a w ~ 
(2 log* + logw B )/(logr x t 2 ) for t < iwhitc and = 
after t w hite- 

We now concentrate on the distribution of light gray 
values corresponding to white strips that have just been 
compressed below the cut-off scale We propose to 
approximate the measured value C as the average of the 
biggest white strip with width A < % , and mixed " gray" 
fluid whose concentration is close to the most probable 
concentration C g . A box with a white strip of scale A thus 
bears a concentration C\ = C g (l — X/wb) and we can re- 
late P{C) to the distribution of widths of the images of 
the injected white strips Q(A). A white strip injected at 
to is transformed into 2 t_t ° images with scales dt o r*~ to . 
Therefore Q(A) = (A/d to ) lo s( 2 )/ lo s( r ) x (1/AlogT), re- 
sulting in 
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P(C) thus has a power-law tail in the light gray lev- 
els whose exponent depends on the mean stretching T. 
We observe satisfactory agreement between this predic- 
tion and both experimental data and numerical 1-D sim- 
ulations (see Fig. Q] (c) and Fig. [3]) where for this tail 
P(C) cx (C- Cg)~ Q with a < 2, consistent with T < 0.5, 
a rather homogeneous stretching. Also note that the am- 
plitude of the light-gray tail decreases as a power law 
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where a(T) = 1 — log 2/ log T, and C m - ln is the small- 
est concentration observed (C m ; n = for t < t w hite and 



C mill = C g (l - d t /w B ) for t > twhite)- For t < t v 



hite 



the 



integral is constant and <r 2 G cc g(t) (x d t cc t~ 2 . On the 



other hand, for t > t 



white 



git) 



2 + a(r) 



For a(r) ~ 2 as we observed, the exponent in the above 
power law is about —4. 

Finally, the dark gray tail corresponds to boxes con- 
taining images of mixed regions from early times - thus 
with an important percentage of black - that have ex- 
perienced little stretching. It is therefore not sufficient 
to consider only the mean stretching T as before, since 
stretching histories far from the mean are involved. Us- 
ing the large-deviation function S for the finite-time Lya- 
punov exponents distribution [l4j |. we derive 



P(C,t) = 



2*-*° exp[-(i-t )S(- 



log W E 

t-to 
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where to(C) is the earlier time at which the mixed region 
had a mean concentration C, so that C = (C)(l — d to ) _1 . 
At a fixed time the dark gray tail decreases as the domi- 
nant contribution (C— (C)) -2 , however the probability to 
observe a fixed concentration value decays exponentially, 
allowing us to neglect a 2 DG in ([T]). Both the (C — (C))~ 2 
shape and the rapid fall-off of the dark-gray tail can be 
observed on Fig. [T] (c) and Fig. [3l 

We now sum these contributions to obtain er 2 (C). 
In the experiment, the crossover £ w hite is estimated as 
30 periods. However, 3-D effects inside the fluid pre- 
vented us from conducting experiments for more than 
35 periods. For this early regime, fitting the data with 
a^y oc (21ogi + logWB)A 2 (black line on Fig.[T](b)) gives 
good results, except close to £ w hite where the contribu- 
tion of the light gray tail starts to dominate. In contrast, 
in numerical simulations we observe (Fig.Q](b)) both the 
(21og£ + logwe)/* 2 behavior (black line) , which can be 
interpreted as in the experiment, and the t~ 4 decay after 
^whitc (100 periods for the case studied) given by a\ G . 

Note finally that the algebraic nature of d{t) permits 
crude "first-order" approximations such as considering 
only the mean stretching given by the Lyapunov expo- 
nent as we did. The injection process dominates other 
mechanisms put forward to analyze concentration dis- 
tributions, such as the evolution of the distribution by 
self-convolution due to the random addition of concen- 
tration levels lH, 16} . The strange eigenmode formalism 
also fails to describe this nonasymptotic regime as the 
spatial mixing pattern is still evolving. 

In conclusion, we propose a scenario for 2-D mixers 
with a chaotic region that extends to fixed walls. As 
soon as the flow is sufficiently slow close to the wall - a 



rather generic situation - and keeps for long times an un- 
mixed pool close to a parabolic point, the unstable man- 
ifold of which feeds fluid into the mixing region, the en- 
tire concentration field is affected, regardless of distance 
from the walls. The algebraic scalings for the variance 
and concentration distributions can be predicted from 
the filamentary stirring pattern generated by a combina- 
tion of stretching, folding, and injection of fluid from the 
fixed walls. Our reasoning could be extended to other 
algebraic expressions of d{t) resulting from different hy- 
drodynamics at the wall. In the present case (no-slip 
wall), we derive a power-law (asymptotically t~ 4 ) evo- 
lution for the decay of the concentration variance, and 
find very good agreement between our analytical, exper- 
imental, and numerical results. Our simulations of the 
well-studied viscous blinking vortex flow Q also yield the 
same algebraic decay for the variance of a coarse-grained 
concentration, obtained by advecting a large number of 
points. These numerical results will be presented in a 
future paper, where we will generalize our approach to 
other systems characterized by continuous injection of 
inhomogeneity, such as open flows. 
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